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By adopting a concept from signal processing, instead of starting from the correlation 
functions which are even, one considers the causal correlation functions whose Fourier 
transforms become complex. Their real and imaginary parts multiplied by 2 are the 
Fourier transforms of the original correlations and the subsequent Hilbert transforms, 
respectively. Thus, by taking this step one can complete the two previously needed 
transforms. However, to obviate performing the Cauchy principal integrations required 
in the Hilbert transforms is the greatest advantage. Meanwhile, because the causal 
correlations are well-bounded within the time domain and band limited in the 
frequency domain, one can replace their Fourier transforms by the discrete Fourier 
transforms and the latter can be carried out with the FFT algorithm. This replacement is 
justified by sampling theory because the Fourier transforms can be derived from the 
discrete Fourier transforms with the Nyquis rate without any distortions. We apply this 
method in calculating pressure induced shifts of H 2 0 lines and obtain more reliable 
values. By comparing the calculated shifts with those in HITRAN 2008 and by screening 
both of them with the pair identity and the smooth variation rules, one can conclude 
many of shift values in HITRAN are not correct. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

In order to carry out forward modeling of atmospheric 
radiative transfer processes and to obtain information on 
the abundances of molecular species, temperature-pressure 
profiles, and other atmospheric properties, one needs accu- 
rate spectroscopic data. This includes not only line positions, 
strengths, and the temperature-dependent Lorentzian half- 
widths, but also pressure induced shifts. Because the 
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ambient atmospheric species, temperatures, and pressures 
are not always amenable to laboratory measurements, or 
because of the large number of transitions possible, one 
often has to rely on theoretical calculations. This is espe- 
cially true for the induced shift because their values are 
smaller than values of the half-width so that to perform 
reliable laboratory measurements becomes more difficult. 

It is well known that there are formalisms available for 
many years such as the Anderson-Tsao-Curnutte (ATC) 
theory [1,2], the Robert-Bonamy (RB) theory [3,4] and 
others with which one can calculate half-widths and 
shifts of molecular spectral lines. These formalisms are 
based on two basic approximations: the binary collision 
and the impact approximations. With these two approx- 
imations, calculations of the pressure broadened half- 
widths and induced shifts of the absorber molecular lines 
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are reduced to thermal averages of the “connected” opera- 
tor Mc(<m) introduced by Fano over all possible collision 
processes involving a pair of molecules [5], and in compar- 
ison with the time of interest, these collisions are assumed 
to be completed instantaneously. In addition, both the ATC 
and RB formalism are semi-classical theories. In other 
words, the translational motion of the pair is treated 
classically, while their internal degrees of freedom are 
treated quantum mechanically. This semi-classical method 
is valid in calculating molecular spectral lines for tempera- 
tures of interest in atmospheric applications. However, the 
current theories are also based on other assumptions. 
Some of these lack sound justification and could lead to 
uncertainties that affect the reliability of the calculated 
results. Besides, it is commonly accepted by most research- 
ers that in comparison with calculated half-widths, calcu- 
lated shifts usually contain larger uncertainties. According 
to our knowledge, there is lack of systematic studies to 
explain why the theoretically calculated shifts are poorer 
than the calculated half-widths. One possible reason could 
be severe cancelations of contributions to the lineshifts 
occur in calculation processes. To identify other reasons 
requires individually investigating the entire calculation 
processes within the specified formalisms. 

With the formalism developed from the introduction 
of the coordinate representation in Hilbert space, [6,7] the 
usual route in calculating converged half-widths and 
lineshifts consists of several steps. For simplicity, we omit 
subscripts attached to following functions. First, one 
calculates the irreducible correlation functions F(t) of 
the S matrix based on the potential and trajectory models 
selected. Then, by carrying out the Fourier transforms of 
F(t), one obtains functions H(co). Next, by evaluating the 
Cauchy principal integrations for H(a>), one obtains l(co), 
the Hilbert transforms of H(co). After these H(co) and 1(a ) ) 
are available, one calculates the real part of S 2 , the 
second-order of the perturbation expansion of the S 
matrix, from H(co) and the imaginary part of S 2 from 
1(a)) for specified lines. Finally, one derives their calcu- 
lated half-widths and shifts. 

Following this route, the entire process goes smoothly 
until one has to do the Cauchy principal integrations. 
There are subroutines available in computer libraries with 
which one can carry out the Cauchy principal integra- 
tions. In general, these subroutines work well for most 
cases. However, it turns out their unstable performances 
do happen occasionally. Given the fact that in practical 
calculations one needs to evaluate a huge number of the 
Cauchy principal integrations of the order of 10 7 , even 
seldom occurring performance failures could cause severe 
consequences because it is not easy to identify all of them 
and to make remedies accordingly. As a result, our 
previous attempts to calculate converged results for line- 
shift stalled at this stage. In fact, except for reporting 
calculated half-widths of H 2 0 lines [7], we have not 
published any calculated shifts for H 2 0 lines because of 
the instability problem mentioned above. Even for our 
calculated half-widths, in order to avoid involving 1(a)) we 
had to follow other researchers’ assumption that the 
contributions from the imaginary part of S 2 can be 
ignored. Fortunately, thanks to the present study the 


validity of this assumption has been verified. Therefore, 
within our formalism it is the Cauchy principal integra- 
tion that causes more difficulty in calculating the lineshift 
than the half-width. Unless we are able to find an 
alternative route to reach /(to), we simply cannot provide 
reliable lineshift values at all. Despite realizing this 
problem for years, we had not found an effective solution 
to overcome this challenge until now. 

Very recently, an idea occurred that because the Four- 
ier and Hilbert transforms are commonly used in many 
science and engineering fields, why not pursue a search for 
tools used in other disciplines? This effort paid off quickly, 
and we found a useful tool that we were looking for. By 
adopting the concept of the causal function from signal 
processing, we have found that instead of starting from 
the usual correlations which are even functions and whose 
Fourier transforms are real, we can consider the causal 
correlations which are neither even nor odd and whose 
Fourier transforms become complex. But the most striking 
feature of the latter’s Fourier transforms is that their real 
and imaginary parts are simply equal to H(m)/2 and I(a>)/ 2, 
respectively. With this tool, the Cauchy principal integra- 
tions are completely obviated. Meanwhile, thanks to sam- 
pling theory, the Fourier transforms can be replaced by the 
discrete Fourier transforms with the Nyquist rate without 
any distortions [8]. Furthermore, the latter can be easily 
carried out with the fast Fourier transforms (FFT) [9]. Our 
numerical tests demonstrate accuracies achieved for H(co) 
and I(co) derived are well satisfied. Especially for 1(a)), 
where one is able to achieve very high accuracies in large 
ranges of a). Thus, with this simple method, the challenge 
has been completely overcome. 

The present manuscript is organized in the following 
way. In Section 2, we only outline the basic of the 
modified RB formalism developed in the coordinate 
representation [6,7], but we provide explanations in detail 
about the correlation functions, their Fourier transforms, 
and the subsequent Hilbert transforms. Section III is 
devoted to the new route in which we start by introdu- 
cing the concept of the causal function and its Fourier 
transform. Then, we briefly explain the sampling theory 
with which one is able to replace a Fourier transform by a 
discrete Fourier transform. Finally, by presenting exam- 
ples, we demonstrate how the new method works in 
evaluating the Hilbert transform without performing the 
Cauchy principal integrations and explicitly show how 
high accuracies of the new method can be achieved. In 
Section 4, we present our calculated shifts for H 2 0 lines in 
the pure rotational band and make comparisons with 
those listed in HITRAN 2008, and other fruitful results 
from the present study are also presented. The last 
Section 5 consists of a brief discussions and conclusions. 

2. General formalism in calculating N 2 induced shifts for 
F1 2 0 lines 

2. /. The RB formalism applicable for complicated potential 
models 

Among all line shape theories available at present, the 
RB formalism is one of the most widely adopted. In fact, it 
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is this formalism with which theoretical calculations of 
the half-widths and shifts of H 2 0 lines have been carried 
out for years [10-13] and the results have played impor- 
tant roles in updating the HITRAN database [14]. We 
briefly outline the formalism used to calculate the half- 
widths and shifts of H 2 0 lines. 

The main computational task for calculating the Lor- 
entzian half-widths and shifts is the evaluations of matrix 
elements appearing in the perturbation expansion of the 
Liouville operator S ( =S, • Sj, where S, and S F are scattering 
matrices in Hilbert space). Usually, in practice, these 
evaluations are limited to the second-orders of the expan- 
sion. Within the RB formalism, the original expressions 
for the pressure broadened half-width and induced shift 
are given by 

7rb = 2ncj 0 Vf(y)dV J 0 2nbdb 

x < l-cos[Si(b) + IrnS 2 (b)]e~ ReS2(f,) >j 2 , (1) 

and 

„ r + OC r+OO 

<5rb = ^ J vf(v)dv I 2nbdb 

x <sin[S 1 (b) + ImS 2 (b)]e~ ReS2(b) > J - 2 . (2) 

where n b is the number density of the bath molecule, f[v) 
is the Maxwell-Boltzmann distribution function, S! and S 2 
are matrix elements in the Liouville space associated with 
the first- and second-orders of the perturbation expansion 
of the Liouville operator S, and < •■■ >j 2 means an average 
over the quantum number j 2 of the bath molecule. 
However, we have found a subtle error in the RB formal- 
ism [15]. After remedying this derivation error, the correct 
expression for the half-width and lineshift in the “mod- 
ified” RB formalism (MRB) become [15] 

„ r+OO r+OO 

7mrb =2^ Jo v f (v)dv J 0 2nbdb 

X { 1 -cos[ < S, (b) > j2 + 1m < S 2 (b) >j2 ]e- Re <^® >i } . (3) 

and 

r +oo 

v/(v)dv / 2nbdb 
Jo 

X Sin( < S, (b) > j2 + Im < S 2 (b) > j2 )e- Re < s ^ 6 »« . (4) 

We note that the essential difference between these 
expressions is that in Eqs. (1 ) and (2) the summation over 
j 2 is outside of the cumulant expansion while in contrast, 
in Eqs. (3) and (4) it is inside. It is worthwhile emphasiz- 
ing that because the bath average is carried out in the line 
space, an average over j 2 must be included. Thus, by 
understanding that Si(b) and S 2 (b) are the bath averages 
associated with the first- and second-order expansions of 
the Liouville operator S, one can suppress all of these 
< ■•• >j 2 in Eqs. (3) and (4). For simplicity, we will omit the 
subscript MRB of y and 8. 

Usually, people prefer to represent S 2 by three compo- 
nents labeled by S 2iOUter /r c ), S 2 , outer /r c ), and S 2 , mldd , e (r c ), 
respectively. In the present study, we follow the same 
custom. As an example, an expression for S 2ou[eri (r c ) is 



given by [6,7] 


5 2,ouCer,i(Tc ) 

1 


X f°dt f dte Ka W' +a ’* xt - ti 
J — oo J— oo 


X 1 V( R (t)) 

x 1 V CR (O) L/iT { m J 2 m 2 ^ , 


( 5 ) 


where R ( t ) represents the translational motion of two 
interacting molecules. In the above equation, p Jz is the 
density matrix of the bath molecule, (m) means summa- 
tions over all magnetic quantum numbers, coy = 
[E (h) (j 2 )-E (b) (j' 2 )]/h, and \j t Timj 2 m 2 > = \j i 'tim i > ® |j 2 m 2 > 
where IjiTiirii > and \j 2 m 2 > are the basis of Hilbert space 
for the H 2 0 and N 2 molecules, respectively, and t is a 
simple notation for k a -k c in |jTm > . In practical calcula- 
tions, how to evaluate these three S 2 terms more accu- 
rately presents a big challenge. 

In order to evaluate the potential matrix elements in 
Eq. (5), one prefers to express V(R(t )) in terms of the 
standard spherical expansions 


V(R(t))= U(L 1 L 2 L;K 1 ;R(t )) 

Li K\L 2 L 

X ]T C(LjL 2 L,m i m 2 m)D^* Ki (Q a )D^* 0 (Q b )Yf m (o}(t)), 

m i mzm 

( 6 ) 

where R (t) is described by its polar coordinates R(t) and 
co(t), and Q a and Q b are orientations of the absorber and 
bath molecules, respectively. In the present study, we 
adopt potential models for the H 2 0 and N 2 molecules 
consisting of long-range dipole-quadrupole and quadrupole- 
quadrupole interactions V dq (Q a ,Q b , R) + V qq (Q^,Q b ,R) 
and a short-range interaction V atom _ atom (£2 d ,£2(,, R ) mod- 
eled by the site-site model. For V dq and V qq , their spherical 
expansions are well known and numbers of their expan- 
sion terms are very limited. For V dq , the expansion con- 
tains only one term (i.e., Lj = l, K i = 0, I 2 = 2, and L=3) 
which varies with R as R~ 4 . Meanwhile, for V qq it has three 
terms (i.e., L i = 2, /Ti = 0, ±2, L 2 = 2, and L=4) which vary 
as R 5 . 

With respect to V a[om _ atom , the expression for the site- 
site model is given by 

( a V (T 6'1 

Vato m-atom&a.Qb. R ) = ^ ~ \ - ( 7 > 

iea jeb l T ij T ij J 


where and By are parameters and ry are distances 
between the i-th atom of the absorber molecule a (i.e., 
H 2 0) and the j-th atom of the bath molecule b (i.e., N 2 ). In 
order to evaluate its matrix elements, one needs to 
rewrite V atom _ a tom in terms of the standard spherical 
expansions 


^atom-atomi^ar^b' ^ (0) 

_ bJ(L\K\L2L,Tlyj},WQ) 

~ 2-^t nL ! +L 2 +q + 2w ft s 

LjKjL 2 Ln {ij) wq K K L ) 


X C(LiL 2 L,tn 1 m 2 m)D%* Ki (G a )D%* 0 (Q b )Yi m (w(m, 


( 8 ) 
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where n {I j} runs over all pairs of atoms in Eq. (7), q = 6 or 
12, w is an integer index from 0 to infinity, and the 
definitions for U(LzKjL 2 L,n[y},wq) can be found in the 
literature [16]. 

As shown in Eq. (8), when one evaluates the matrix 
elements of V a[om - atom (t), there would be a lot of summation 
indices. In order to evaluate S 2outerl {r c ), one needs to con- 
sider the products of two potential matrix elements. Thus, 
the total number of summation terms in Eq. (5) becomes so 
large that the evaluations are intractable unless cut-offs are 
enforced. It turns out that there are two kinds of cut-offs 
associated with the summation indices in Eqs. (8) [6,7]. The 
first is a cut-off to limit the set of irreducible tensor ranks 
and L 2 and the second is to set an upper limit for the index 
w. In the literature, if one chooses 8 as the maximum of 2w, 
the cut-off is said to be the 8-th order. It is the introduction 
of these two cut-offs that opens the possibility that the 
results derived are not converged. We have provided in 
detail analyses on the convergence problems existing in this 
method and we do not repeat them here. Readers interested 
in this subject are referred to our previous works [6,7]. 


2.2. Irreducible correlation functions of the S matrix in the 
coordinate representation 


In order to overcome the convergence problem, we 
introduce the coordinate representation [6] whose basic 
idea is briefly outlined here. It is well known from quantum 
mechanics that one has the freedom to choose different 
representations in Hilbert space. The state representation 
that is commonly used consists of the basis set ® 

\j 2 m 2 > which are eigenfunctions of the Hamiltonian of the 
molecular pair. On the other hand, one can select the 
orientations of the pair of molecules as the basis set in 
Hilbert space; i.e„ |(5(f2 a -f2 aa ) > ® \d(Q b -Q ba )} where 
Q aal and C2 bal represent orientations of H 2 0 and N 2 specified 
by a, respectively. By choosing the coordinate representa- 
tion, the potential becomes a diagonal operator and can be 
treated as classical functions [6,7]. Then, by introducing the 
correlation functions described later, one is able to calculate 
the half-widths with high accuracy, no matter how com- 
plicated the potential models are and no matter what kind 
of trajectories are chosen. In any case, the convergence of 
the calculated results is always guaranteed. 

We will not provide a detailed development of the new 
formalism here and readers can find its derivations from 
our previous works [6,7]. We only present some important 
formulas here. The correlation functions are defined by 


/ OO 

dt'G t]Kl jq i2 (t , + t/2,t'— 1/2), 

-OO 


(9) 


where G 1]KlJf;i2 (t,t') is given by 

GL,K,K',L 2 (t,t ) - 47lfi 2 (2il+1) 2 (2i2 + 1) 2 

x(— if 1 £ (—1/6 +G +9(21+ 1) 

L 

x u(Lii 2 L; /Ci ; R( t))u(L]L 2 L; I<\ ; R(t'))P L (c os0 [t '). 

( 10 ) 

In Eq. (10), the factor (-l) (t| +L i +L > — l because the summa- 
tion index L must satisfy L 1 +I 2 +L=even, and are 


angles between two vectors R ( t ) and R (t'j. The set of 
indices used to label the correlation functions consist of one 
tensor rank L] with two subsidiary indices If, Kf related to 
H 2 0 and another tensor rank L 2 for N 2 . Because N 2 is a 
diatomic molecule, L 2 must be even. If one chooses the II R 
representation to develop the H 2 0 wave functions where 
two H atoms are symmetrically located in the molecular- 
fixed frame, values of K, and Kf must also be even. 

In the standard method, the functions required to be 
evaluated are the resonance functions and their total number 
depends on how many combination choices of summation 
indices L t , /<j, L 2 , L, and q, n^y w are used in Eq. (8) with the 
restrictions enforced by the two kinds of cut-offs. As these 
two cut-offs increase, the number of resonance functions 
increases very quickly and to evaluate all these resonance 
functions become formidable. On the other hand, with the 
new method, the functions required to be evaluated are the 
correlation functions and their total number is determined 
only by the upper limits of Lj and L 2 restricted by the first 
kind of cut-off [6,7], This means that no matter how high the 
second kind of cut-off is, the number of correlation functions 
remains unchanged. In addition, because V a[om _ a[om becomes 
ordinary functions, one can choose any order cut-off one 
wants. No matter how high this cut-off order goes, there are 
only small differences determined by the needed computa- 
tional resources. As a result, we can use both cut-offs which 
are sufficiently high to guarantee full convergence in practical 
calculations. It is easy to show that if one sets 2 as the 
maxima for both L ^ and L 2 , the number of correlations is 20. 
Meanwhile, if one increases the maximum for U from 2 to 3 
or 4, the number of correlations increases from 20 to 38 or 
88, respectively. Finally, by setting 4 as the maxima for both 
Li and L 2 , the number becomes 132. 

Then, in terms of the correlations, one is able to 
rewrite Eq. (5) as 

POO 

$2, outer, iO"c) = E L dt K*^ t)W '>F L x KjK\L 2 (t)- 

l,k 1 k\l 2 J0 

( 11 ) 

In deriving the above expression, two functions which are 
independent of the potential and trajectory models are 
defined by 

Wf^ K , (t) = J2(2fi + 1 JDO'iTjfo; I, JC, ) 

j'i X i 

xD(j i rJ' i T' i ;L,K\)e“ 0 ^ t , (12) 

where DQ't/t';/ K ) is defined by 

D(jzfT'-LK) m £ (-V k U[ z U i ;_ KT ,C(jj'L,kI<-kl<), (13) 

k 

and 

W[ h >(t) = + 1 X 2 l 2 + 1 )p j2 C 2 (f 2 j' 2 L2,000)e m ‘^, (14) 

J2J2 

respectively. In Eq. (13), lfi ki . are expansion coefficients of 
the H 2 0 wave functions in the II R representation over the 
symmetric top wave functions | jkm} with k=—j, — j+1, 
j such that 

\jxm > =J2 u iz \J km > ■ 

k 


(15) 
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2.3. Fourier transforms of the correlations and subsequent 
Hilbert transforms 

As further steps in developing S 2 outer.i(r c ) from Eq. ( 11 ), 
we apply two tools, i.e., the Fourier transform and the 
Hilbert transform. First of all, by introducing H LjKiK ' jL2 (m), 
the Fourier transforms of the correlation functions, 
defined by 

1 T°° ■ , 

Hl,k,k\l 2 ( ( °) = eltu F 1]K]JC ' ]l2 (t)dt, (I®* 

one is able to obtain an expression for the real part of 

■^2 ,outer,i(tc) as 

^^2,outer,i(tc) = ifw ^ . /* (2j,- + * ) 

V z l,k,k\l 2 jft 

X Lj Kj )D(j i z i f i x' i ; Lj I <\ ) 

x £(2 j 2 + 1 )(2j 2 + 1 )Pj 2 C 2 (j2j 2 L 2 ,000) 

J 2 J 2 

x ^L 1 K,K , t L 2 ( a> j i z i jir' + COj 2 J' 2 ). ( 17 ) 

Secondly, by introducing Il 2 k 2 k\l 2 (0)< the Hilbert trans- 
forms of H 1]K] , C ' l2 (co), defined by 

k,K,K-L 2 m=^pJ^ dof - c ^H LlKlK . L2 (of), (18) 

where P means the principal part, one is able to obtain an 
expression for ImS 2 ou t er ,i('r c ) which is the same as ReS 2-ou _ 
ter,i( r c) except for a replacement of H LlK ^ L2 (m w . t ,+a) j2f2 ) 
by h,K,K\L 2 (o>] ttAti +o) ]2j ' 2 ) in Eq. (17). With respect to 
expressions for ReS 2 , 0 uter/r c ) and ImS 2-out er,X r c). they are 
same as their corresponding partners except for the 
replacement of the initial quantum numbers j, t;, and 

T'j by the final ones jf t f jf, and t f. Meanwhile, the 
expression for S 2 middle (r c ), which is real, is given by 

S 2 .middie(r c ) = s/2n J2 {(-l) 1, (2j, + l)(2j / + l) 
l,k,k\l 2 

x W(jij f jj f ; lLj )D(j i Tj i x i ; L] Kj )D(j f r f j f T f ; Lj Kj )} 
x ^^(2 j 2 + 1 )(2 ,j 2 + 1 )Pj 2 C 2 Oz) 2 i 2 ,000) 

J 2 J 2 

xH L 1 K 1 K' 1 L 2 ( ffl j 2 j' 2 )' (19) 

After all these expressions are available, we have obtained 
two basic formulas, 

ReS 2 (r c ) = ReS 2outeri ( r c ) + ReS 2outl , r ^(r c ) + S 2mi(id ( e ( r c ), (20) 

and 

ImS 2 (r c ) = -ImS 2 , 0Utf , rii (r c ) + ImS 2 , 0Uteri/ (r c ), (21) 

where the minus sign of Im S 2-out er,;(r c ) results from the 
fact that S 2 (r c ) = S 2 ,oiiter,i(r c ) ^ S 2 ,outerj( l~c) "F S 2 ,middl(‘(. h.) [10]. 

It is worth mentioning that the relationship among 
these functions is that H LjKiKjL2 (cl>) are the Fourier trans- 
forms of F llK]K![l2 (t), meanwhile I L;KtK \L 2 (0J) are the Hilbert 
transforms of H L]KiK ' jL2 (oj). In addition, because F I|K|K . t2 (t) 
are even functions, their Fourier transforms H LjKjK - iL2 (o j) 
are real and even, and the subsequent Hilbert transforms 
I/., k, k\ l 2 ( 0) ) are real and odd. Finally, we note that all these 
functions F t]K]ril2 (t), H LlKlK - L2 (o)). and / tlKl r,L 2 (ffl) are 
associated with specified trajectories. As a result, all of 


them depend on r c . For simplicity, we have omitted the 
argument r c in their notations. 

In practice, however, this two steps method may not 
be the best way to follow. For the first step, usually the 
continuous Fourier transforms are replaced by the dis- 
crete Fourier transforms with proper samplings. The latter 
can be easily derived with the FFT algorithm [9] which is 
cited as one of the most significant contributions to 
numerical analysis of the 20-th century [8]. As a result, 
there is no obstacle in carrying out this step. However, the 
second step requires one to evaluate the Cauchy principal 
integration. Although subroutines for carrying out the 
Cauchy principal integration are available in main com- 
puter libraries, their performances are not always satis- 
factory. In general, these subroutines work well for most 
of cases. However, their unstable performances do happen 
occasionally. There are a lot of independent Hilbert trans- 
forms fi,K,iqi 2 (®) and if one considers 132 correlations 
this number is 39. In addition, all these functions have 
two dimensions. In order to depict each of them well, one 
has to use high resolutions by selecting several hundred 
points for each of its two variables such that the corre- 
sponding range is well covered. As a result, in practical 
calculations sometimes one needs to evaluate a huge 
number of the Cauchy principal integrations in the order 
of 10 7 . An obstacle here does not result from the compu- 
tational burdens, but from unstable performances of the 
subroutines occasionally happening because one is not 
able to pick up all unstable results from such huge data 
pool and to make modifications accordingly. As a result, 
the ability to derive reliable two dimensional profiles of 
h.,K,K',L 2 ( Cj) ) becomes a big challenge. In fact, in spite of 
that the formalism in calculating the half-widths and the 
lineshifts based on the coordinate representation has 
been developed several years ago, [6] we have not 
reported any calculated results involving evaluations of 
h,i qk\l 2 ((o). We have presented calculated half-widths for 
the H 2 0-N 2 system [7], but we have assumed there that 
contributions from ImS 2 (r c ) to calculated half-widths can 
be ignored. Meanwhile, we have not reported any calcu- 
lated shifts at all. In summary, in order to make progress, 
to find an alternative way to evaluate / /|K]/f - t2 (oj) becomes 
mandatory. Recently, we began to wonder whether taking 
the two steps is the only way to find Il,k^k\l 2 (. co ') or can 
one derive them directly from F L]K]K]t2 (t). 

3. Causal correlation functions and Fourier transforms 

3.3. Causal correlation and Fourier transform 

We briefly explain a concept about the causal function 
and related topics which have been widely used in signal 
processing and many other engineering and science dis- 
ciplines. First of all, let us assume F(t) of interest is an 
even function. We want to find its Fourier transform H(cu) 
together with the subsequent Hilbert transform I(a>). The 
usual route which is depicted in Fig. 1 is to cany out the 
Fourier transform of F(t) first, then one performs a sub- 
sequent Hilbert transform of H(at ) to derive /(co). A 
question arises that is there a way to establish a direct 
link between F(t) and I(co )? 
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Fig. 1. A diagram to show the usual route to derive the Fourier transform 
H(co) from the function F(t) and a subsequent Hilbert transform l(c o) of 
H(ca). Is there a way to establish a direct link between F(t) and l(c u)? 


Because the Hilbert transform including the Kramers- 
Kronig relations associated with it has been applied in 
solving many engineering and science problems, we 
believe that by searching tools used in other disciplines, 
one may be able to find a useful one. Fortunately, the 
useful tool does exist in other fields such as signal 
processing [17]. By adopting this method, instead of 
starting from the function F(t ) itself one defines its causal 
function denoted by F(t) as 

F(t) = F(t)0(t), (22) 

where d(t) is the unit step function (the Heaviside step 
function) defined by 
0(t) = 1 t > 0, 

0(t) = 0 t<0. (23) 


Then, one considers the Fourier transform of the causal 
function F(t) which can be expressed by a convolution of 
the Fourier transforms of its products, i.e., H(co) and 0(co), 


1 r°° ~ ■ 1 r°° 

H(cd) = / F(t)e' wC dt = / Ft(CD')0(CD-CD')dcD' . 

V2n J — oo V2n J —oo 

(24) 


It is worth mentioning that because F(t) is neither even 
nor odd, H(co) is a complex function. On the other hand, 
the Fourier transform of 0(t) is well known and can be 
expresses as 


0(cd) = 


1 

V2n 


(nd(co)+iP—), 

CD 


(25) 


where P means to take the Cauchy principal value upon 
integration. With Eq. (25), one can rewrite Eq. (24) as 


H(co ) = 


/ OO 

H(co' 

-OO 


nS(co—CD')+iP - 
H(co') 


dco' 


H(co)+i-P f°° (j C0 ' 

K J - co CO-03' 


(26) 


With Eq. (26), one can conclude that H(ct>) and /(co), the 
Fourier transform of the function F(t) and the subsequent 
Hilbert transform, equal to twice the real and imaginary 
parts of the Fourier transform of the causal function F(t), 

H(co) = 2R eH(co), 

I(co) = 21m H(co). (27) 


As a result, by just taking one step, one is able to derive 
H(m) and I(ca) simultaneously. In comparison with the 


usual method, the new method enables one to reduce two 
steps to one. But the most important advantage of the 
new method is its ability to evaluate Cauchy principal 
integrations without carrying out these integrations 
themselves at all. As mentioned above, thanks for the 
powerful FFT, the remaining step can be well accom- 
plished. In fact, the computer codes used here are the 
same as those used in completing the first step of the 
previous method. 

Of course, the success in applying the new method here 
depends on how accurately the Hilbert transform /(co) can 
be obtained. Before we present test results to demonstrate 
accuracies of this method could achieve, we briefly explain 
a technique used to derive the Fourier transform of a 
function F(t). This technique has been successfully used to 
derive FI(co) in our previous half-width calculations [6,7] 
and will be used in calculating H(co) here. 

3.2. Sampling the correlations and calculating the Fourier 
transforms 

The method used to derive the Fourier transform of a 
function is based on the sampling theory. Readers can find 
detailed explanations about this technique in the book by 
Weaver [8[. The basic idea is by sampling, one converts a 
function of interest to a sequence, calculates its discrete 
Fourier transform, and then relates this transformed 
sequence to the Fourier transform of the original function. 

When one samples a function F(t ) with a sampling rate 
At, one converts F(f) to a sequence (F(n)} whose terms are 
values of F(t) at the discrete locations t=nAt. According to 
the Whittaker-Shannon sampling theorem, if the function 
F(t) is band limited with handwidth Qf and the sampling rate 
is chosen such that At=l/(2/2/) called as the Nyquist rate, 
then F(t) can be recovered uniquely and exactly from the 
sampled sequence (F(n)}. We note that a function is called 
“band limited with bandwidth Qf if its Fourier transform 
H(tu) has bounded support over the interval [ - Qf, Qf], 

Usually both the frequency f and the angular frequency 
co are referred by people as the frequency. In order to 
avoid confusion, we use the notation Qf to indicate the 
one associated with the frequency f, but not the angular 
frequency co. However, because many spectroscopists 
prefer to use the angular frequency co in their studies, 
we follow their preferences here. Then, the interval 
should be understood as \ —2n Qf 2n Qf], 

One assumes the function F(t) has bounded support 
over the time domain [ - T, T] and as well as being band 
limited over [ — 2nQf 2nQf\. Because the correlation 
functions introduced in calculating the half-widths and 
lineshifts are smooth functions and their magnitudes 
decrease very quickly to zero as t increases, we believe 
this assumption is applicable for all of them. First, one 
samples F(f) with the Nyquist rate Af=l/(2/2/) and 
chooses the number of samples such that NAt > 2 T and 
obtains a sequence (F(n)}. According to the sampling 
theory, the original function F(t) can be completely 
recovered from (F(n)} such that 

N-l 

F(f)= F(n)sinc(2nQf(t-nAt)). (28) 

n = 0 
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Then, by taking the Fourier transform both side of the 
above equation, one can obtain 

i N— 1 

H(co) =^ry F(n)e~ icjnAt . (29) 

2U f n = 0 

Next, one samples H(co) with Acu=27t/(2T) (i.e., with the 
Nyquist rate A/=l/(2T)=2f2j/N) and obtains another 
sequence in the frequency domain 

1 N -> flN-l 'I 

H(m) = F(,n)e- imnAmAt = N\t ■ e- 2nimn ' N \ 

2Q f n = 0 l N ^0 J 

= NAt-T[F(n)}, (30) 

where 2 r {F(n)} denotes the discrete Fourier transform of 
(F(n)}. With Eq. (30), one can conclude that the sequence 
{H(m)} are equal to the discrete Fourier transform of the 
sequence (F(n)} multiplied by NAt. Meanwhile, because 
one chooses the Nyquist rate in deriving the sampling 
sequence { H(m )}, it is guaranteed by the sampling theo- 
rem that the Fourier transform H(o>) can be derived from 
{H(m)} without any distortions. Now, we have completed 
a brief explanation about how to effectively and accu- 
rately derive the Fourier transform for the function F(t). 

in our formalism to calculate the half-widths and line- 
shifts, instead of t and co we prefer to use the dimension- 
less arguments z and k as the argument of the correlation 
functions, and the argument of their Fourier transforms 
and the subsequent Hilbert transforms. According to our 
numerical tests, magnitudes of the correlation functions 
could decrease by about 5-7 orders as z changes from 0 to 
12.8. In fact, when z increases from 0 to 12.8, the distance 
between two interacting molecules have increased from 
the closest distance r c to 12.8 x r c . Roughly speaking, the 
latter is at least about 40 A and at this distance interaction 
between the two molecules is almost negligible. Therefore, 
we can consider the correlation functions are limited in the 


range [ -Z, Z] = [- 12.8, 12.8]. In practice, it is better to use 
finer sampling rates Az and larger numbers N such that 
NAz > 2Z. In the Section 3.3, we will present several 
samples to show how this method works. Finally, it is 
worth mentioning that the range of k accessible from this 
method is [ - rc/Az, 7t/Az] and one has to make sure that the 
Az is small enough that this range would be larger than the 
whole k range of interest. 

3.3. The accuracy check of calculated Flilbert transform 

As the first example, we consider the function F(t) to 
be a Gaussian 

F(t) = e~ (t2 / 2 \ (31) 

because its Fourier transform Ff(to) is also a Gaussian 

H(co) = e~ ( “ 2 / Z) , (32) 

and the subsequent Hilbert transform /(co) is well known. 
In fact, the Hilbert transform of a Gaussian is a function 
called as the Dawson’s integral and it can also be related to 
the imaginary error function. In applying the new method, 
one can introduce the corresponding causal function 

F(t) = e- (t2 / 2) 0(t). (33) 

Then, after carrying out the Fourier transform of F(t) with 
the technique described above, the Hilbert transform of 
H(co) can be obtained from the imaginary part of this 
Fourier transform. Because values of the Dawson’s inte- 
gral are known, by comparing calculated values from F(t) 
with them, one can check how accurate the calculated 
values are. 

By setting T=12.8 and NAf=8T, we have calculated 
the Fourier transform of F(t) and obtained I(co ) in the co 
range of [0,1000] with several different choices of N and 
we list some of the results in Table 1 together with the 


Table 1 

Values of /(to) derived from using different sampling number N. 


CO 

Dawson’s integral 

Calculated /(co) 



N=2097152 
At=0.1 x2~ u 

N=131072 
At=0.1 x2~ 7 

N= 65536 
At=0.1 x 2 -6 

1.0 

0.578290E + 00 

0.578290E+00 

0.578290E + 00 

0.578290E + 00 

10.0 

0.8061 16E- 01 

0.8061 16E- 01 

0.8061 18E — 01 

0.8061 24E- 01 

20.0 

0.399947E — 01 

0.399947E — 01 

0.399951 E- 01 

0.399963E — 01 

30.0 

0.266258E — 01 

0.266258E — 01 

0.266264E — 01 

0.266282E-01 

40.0 

0.199596E — 01 

0.199596E — 01 

0.199604E — 01 

0.199628E-01 

50.0 

0.159641E — 01 

0.159640E-01 

0.159651E — 01 

0.159681E-01 

60.0 

0.133018E — 01 

0.133017E — 01 

0.133029E — 01 

0.133066E — 01 

70.0 

0.114007E — 01 

0.1 14006E — 01 

0.114020E — 01 

0.114063E — 01 

80.0 

0.99751 2E- 02 

0.997502E — 02 

0.997664E — 02 

0.9981 52E- 01 

90.0 

0.886648E — 02 

0.886639E — 02 

0.886821 E- 02 

0.887369E — 02 

100.0 

0.797964E — 02 

0.797954E — 02 

0.7981 56E- 02 

0.798765E-02 

200.0 

0.398952E — 02 

0.398940E-02 

0.399344E — 02 

0.400566E — 02 

300.0 

0.265964E — 02 

0.265964E — 02 

0.266571 E- 02 

0.26841 2E- 02 

400.0 

0.199472E — 02 

0.199475E — 02 

0.200285E — 02 

0.202755E — 02 

500.0 

0.159578E — 02 

0.159581E — 02 

0.160596E — 02 

0.163708E — 02 

600.0 

0.132981E — 02 

0.132985E — 02 

0.134206E — 02 

0.137978E — 02 

700.0 

0.113984E — 02 

0.113989E — 02 

0.1 15416E — 02 

0.119870E — 02 

800.0 

0.997357E — 03 

0.99741 9E- 03 

0.101378E — 02 

0.106537E — 02 

900.0 

0.886539E — 03 

0.88661 0E- 03 

0.905067E — 03 

0.964032E — 03 

1000.0 

0.797885E — 03 

0.797964E — 03 

0.818543E — 03 

0.885226E — 03 
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values of the Dawson’s integral. As shown in the table, the 
method works excellently. In general, the smaller to is, the 
higher the accuracy of I(co). In addition, the larger the N is, 
the higher the accuracy. For the highest choice of 
N=2097152 ( = 2 21 ) and At=0.1 x 2~ n , the error is 
always less than 0.01% in the whole to range. With the 
choice of N=262144 ( = 2 18 ) and At=0.1x2~ 8 whose 
results are not listed here, the errors at ct>=10, 100, and 
1000 are 0.0002%, 0.005%, and 0.64%, respectively. With 
the moderate choice of N= 131072 ( = 2 17 ) and 

At=0.1 x2~ 7 , the errors at m=10, 100, and 1000 are 
0.0002%, 0.024%, and 2.59%, respectively. Finally, using 
the lowest one of N = 65536 ( = 2 16 ) and Af=0.1 x2~ 6 , 
these errors become 0.001%, 0.1% and 11%. With respect 
to H(co), our tests show that the calculated results can also 
match the exact Gaussian function of exp(-co 2 /2) very 
well. More specifically, as magnitudes of H(co) reduce by 5 
orders from its maximum (i.e„ 1 at m = 0), the errors still 
remain less than 1%. 

As a second test, we consider a more realistic problem. 
There are many early studies of broadened half-widths 
and induced shifts for two interacting molecules [16,20]. 
In these works, by considering simple potential models 
consisting of only electric multipole interactions and 
choosing simple trajectory models, such as the straight 
line or the parabolic trajectory, people have introduced 
resonance functions and developed methods to evaluate 
these functions. Thanks to their contributions, those 
resonance functions associated with the dipole-dipole, 
the dipole-quadrupole, and the quadrupole-quadrupole 
interactions are available in the literature [18-22]. With 
our formalism, we consider a pair of two linear molecules, 
choose the parabolic trajectory model and a simple 
potential model consisting of only the dipole-dipole, the 
dipole-quadrupole, and the quadrupole-quadrupole 
interactions. In this case, the correlation functions have 
simple properties. There is no interweaving between their 
distance dependence and their internal degree depen- 
dence. Besides, the multipole-moments can be simply 
factored out from their expressions. By exploiting these 
properties, one is able to introduce simplified correlation 
functions F n (z), F 12 (z), and F 22 (z) which are common for 
all trajectories and are independent of the multipole 
moments. Then, with the new method, by introducing 
the corresponding causal partners Fn(z), Fi 2 (z), and 
F 22 (z), one can find out their Fourier transforms and 
obtain Hu(k), Iu(k), Hi 2 (k), /i 2 (k), H 22 (k), and / 22 (/c) 
accordingly. It turns out that the latter should match the 
corresponding resonance functions exactly. Therefore, we 
can use these resonance functions as benchmarks to 
check how accurate our results are. 

Starting from the causal functions of Fn(z), Fi 2 (z), and 
F 22 (z), we have calculated their Fourier transforms with 
the choice of N = 262144 and At=0.1 x 2“ 8 and the latter’s 
real and imaginary parts are our calculated Hn(k), In(k), 
H 12 (k), Ii 2 (k), H 22 (l<), and I 22 (k). Based on the fact that 
H Lj i 2 (k) are even functions of k and Ii t i 2 (k) ate odd ones, 
one only needs to plot half of them along the positive k 
axis. Meanwhile, because magnitudes of h,L 2 (k ) decrease 
more slowly as k increase than magnitudes of Ff t]l2 (k) do, 
in order to present their profiles more completely one has 


to plot h,L 2 (k ) in a larger range of k. We present calculated 
Ffn(k), Fh 2 (fc), and H 22 (k) together with the corresponding 
resonance functions in Fig. 2(a), (c), and (e). Meanwhile, 
we present calculated / n (k), h 2 {k), and I 22 (k) together 
with the corresponding resonance functions in Fig. 2(b), 
(d), and (f) where the k axis is given on a logarithmic 
scale. As shown in Fig. 2(a)-(f), these H LjL2 (k ) and /,, t2 (/<) 
derived from the causal correlations F 1]l2 (z) match the 
corresponding resonance functions very well. Finally, 
there is a small structure of the resonance function 
associated with the quadrupole-quadrupole interaction 
appearing at around k= 7.5 in Fig. 2(f). This minor 
structure indicates the formula used to represent this 
resonance function contains a small error. 

3.4. Samples ofH L ^ Ll (k) and I L , KlK \L 2 (k) derived from the 
causal correlations F LlK]K - , 2 (2) 

First of all, we would like to emphasize that when the 
potential models contain the atom-atom component, the 
distance and the internal degree dependences of the 
correlations become interweaving. As a result, the corre- 
lations F 1]K]J( ' l2 (z) and their associates depend on the 
parameter r c used in labeling collision trajectories. For 
simplicity, we have often omitted this argument r c in 
expressions for these functions, but we will add it in here. 
Base on the potential model used in updating HITRAN 
2008 [14], choosing the 20-th cut-off, including 132 
correlations, and selecting the “exact” trajectory model, 
we have derived functions of H 1]Kl ^ t2 (k,r c ) and 
h 2 K 2 K' 2 L 2 (k,r c ) from all the causal irreducible correlations 
Fl 2 k,k\l 2 (2Jc) with Li <4 and L 2 <4 which are indepen- 
dent. Because magnitudes of h 2 KiK' t L 2 (k,r c ) decrease more 
slowly as /<-> oo than those of Ff tlK]/ q i2 (k,r c ) do, one needs 
to use larger scales of k in plotting them. 

In order to shorten the manuscript, we only present 
Ff 10 o 2 (k.r c ) and IiooiU^rf) derived from the most important 
causal correlation F W02 (z,r c ) associated with the leading 
dipole-quadrupole interaction of the H 2 0 and N 2 pair and 
other terms with the same symmetry in the spherical 
expansion of the atom-atom interaction in Figs. 3 and 4. 
Because H wo2 (k,r c ) is an even function of k and Iioo 2 (k,r c ) 
is an odd one, one only needs to plot a half part for each of 
them. More specifically, the range of k used in plotting 
hoo 2 (k,r c ) in Fig. 4 is [0,120] which is 5 times the range of 
k (i.e. [0,24]) used in plotting Fli 002 (k,r c ) in Fig. 3. With 
respect to the range of r c , because the minimum of the 
closest distance r cmin derived from the “exact” trajectory 
model at 296 I< is 3.5224 A and the dipole-quadrupole 
potential component represents a long-range interaction, 
we have chosen the plotting range of r c is [3.525 A, 
5.025 A]. By adopting these range choices, the most 
important parts of Hi 0 o 2 (k,r c ) and / 1002 (k,r c ) are covered 
in plotting. 

As shown in Figs. 3 and 4, profiles of Hi 0 o 2 (k,r c ) and 
I 1002 (k,r c ) in the region of r c >4.0A somehow resemble 
those provided in Fig. 2(c) and (d). This implies that in 
this region they are mainly determined by the dipole- 
quadrupole interaction and their magnitudes would 
vary as rf as r c -» 00. On the other hand, in the region 
of r c < 4.0 A their profiles contain more structures. 
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Fig. 2. Calculated Hu(k), /, ,(/;), H l2 (k), U 2 (.k), and H 22 (k), I 22 W from the causal correlations Fn(z),Fi 2 (z), and F 22 (z). They are presented in (a)-(f) by red 
dotted curves. Meanwhile, the corresponding resonance functions available in literary are given by black solid lines. 


These structures represent combinations of contribu- 
tions, not only from the dipole-quadrupole interaction 
and the related terms in the spherical expansion of 
the short-range atom-atom interaction, but also cou- 
plings between them. In comparison with H 1002 (k,r c ) 
shown in Fig. 3 , magnitudes of /1002 (k,r c ) decrease more 
slowly such that one have to use a larger scale in 
plotting Fig. 4 . 


4. Calculated N 2 induced shifts of H 2 0 lines in the pure 
rotational band 

In the present study, we consider H 2 0 lines in the pure 
rotational band. By making an approximation to replace 
the integration over the velocity by the averaged velocity 
v ( = y/8 kT/nm) and changing the integration variable 
impact parameter b to the closest distance r c of the 
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4.275 


4.525 


4.775 


5.025 24 


Fig. 3. Fourier transform H W02 (k,r c ) (in ps -2 ) at T=296 K for a molecular pair of H 2 0-N 2 as a two dimensional function of k (dimensionless) and r c (in A). 
As r c ->cn, magnitudes of H 10 o 2 (k,r c ) vary with r c as r c . This function is derived from the real part of the Fourier transform of the causal correlation 
F 1002 (z,r c ) with the choice of N=262144 and Az=0.1 x 2 -8 . 



Fig. 4. Same as Fig. 3 except for the Hilbert transform Iwo 2 (k,r c ) derived from the imaginary part of the Fourier transform of the causal correlation 
£ 1002 (z,r c ). As r c -» oo, magnitudes of I\ooo(k,r c ) very with r c as r c . 


trajectories, one can obtain simplified expressions for the 
half-width and the lineshift 

y = ^ 27 zb [1 -cos(ImS 2 (r c ))e _ReS2(rc) ] dr c , 

( 34 ) 


and 

S=M-[ 27rb(^-)sin[ImS 2 (r c )] e -ReS2(rc) dr c , (35) 

Z7TC Jr cm [ n Of c 

respectively, where r cmin is the minimum value of r c 
corresponding to strictly head-on collisions. After R eS 2 (r c ) 
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Fig. 5. Comparisons between the calculated N 2 -broadened half-widths 
obtained from excluding and including contributions from lmS 2 (r c ). They 
are plotted by A and x , respectively. The 1639 lines in the H 2 0 pure 
rotational band are arranged according to the ascending order of the 
calculated half-width values without contributions from ImS 2 (r c ). 

and ImS 2 (r c ) of lines of interest are available, one can 
easily derive their half-widths and shifts from the above 
two expressions. 

4.1. Contributions from ImS 2 (b) to calculated half-widths 

Usually, people assume that effects from the ImS 2 (r c ) 
term can be ignored in calculating the half-width and 
Eq. (34) can be simplified as [3] 

y= ITc £>(£) [1 - e " ReS2(rC,]drc ' (36) 

In order to justify this assumption, we have calculated the 
N 2 -broadened half-widths for all 1639 H 2 0 lines in the 
pure rotational band listed in HITRAN from both 
Eqs. (34) and (36) and made comparisons in Fig. 5. Both 
these two calculations are made based on the potential 
model used in updating HITRAN 2008 [14] and with the 
20-th order cut-off, including 132 correlations, and adopt- 
ing the parabolic trajectory model. As shown in the figure, 
the differences between these two results are very small. 
More specifically, among all the 1639 lines there are only 
12 lines with their relative errors above 3% and the 
maximum one (happening for the transition of 17 710 <- 
17 4i13 at 410.205 cm -1 ) is 5.1%. Given the fact that the 
formalisms (i.e., the original and modified RB formalisms) 
themselves contain other uncertainties larger than the 
differences reported above and there is an obstacle to 
accurately derive the ImS 2 (r c ) term, to ignore contribu- 
tions from ImS 2 (b) to calculated half-widths is an accep- 
table and justified approximation. Of course, if one knows 
how to accurately evaluate ImS 2 (r c ) which are necessary 
for calculations of the lineshifts, it is better to take into 
account of contributions from ImS 2 (r c ) to the half-widths 
using the Eq. (34). 

4.2. Calculated induced shifts 

Based on the potential model used in updating HITRAN 
2008 and with the 20-th order cut-off, including 132 
correlations, and adopting the parabolic trajectory model, 


Fig. 6. A comparison between the induced shifts listed in HITRAN 2008 
and our calculated values. They are plotted by A and x , respectively. 
The 1639 lines in the H 2 0 pure rotational band are arranged according to 
the ascending order of the calculated shift values. 



Number of lines 

Fig. 7. Comparisons between the calculated N 2 induced shifts obtained 
from the original RB formalism and from the modified version. They are 
plotted by A and x, respectively. The 1639 lines in the H 2 0 pure 
rotational band are arranged according to the ascending order of the 
calculated shift values with the modified RB formalism. 

we have calculated N 2 induced shifts for the H 2 0 lines in 
the pure rotational band. We present our results in Fig. 6 
together with those listed in HITRAN 2008. As shown in 
the figure, there are significant differences between ours 
and that listed in HITRAN. More specifically, among the 
1639 lines there are 649 lines with the relative differences 
above 50%, 746 lines with the differences between 10% 
and 50%, and only 244 lines with less than 10%. It is worth 
mentioning that most of the values in HITRAN 2008 come 
from theoretical calculations [14]. It implies that these 
two theoretical calculations with the same potential 
model differ markedly from each other. 

We would like to emphasize that in order to make 
comparisons more meaningful, we assumed the same 
potential model used in updating HITRAN 2008. Because 
this potential model is poor [6], we do not think our 
current calculated shift values are good. Therefore, read- 
ers should not consider the predicted results in Fig. 6 as 
reliable “true” values. 

At this stage, we would like to take this opportunity to 
demonstrate the importance to use the modified RB 
formalism instead of the original RB one. We present 
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comparisons between calculated N 2 induced shifts of the 
1639 H 2 0 lines derived from the RB and the modified RB 
formalisms in Fig. 7. As shown in the figure, differences 
between them are pretty large. There are 384 lines with 
relative errors above 30%, 767 lines with 5-30%, and 488 
lines with less than 5%. Because both these calculated 
shift values could be equal to or closer to zero, in order to 
avoid artificially enlarging error values we define the 
relative error here by choosing the one with larger 
magnitude as the denominator. Thus, one can conclude 
that the effects from the modification on the calculated 
shifts are significant. 

4.3. Applying the two rules to calculated lineshifts 

Recently, based on the properties of the energy levels 
and wave functions of H 2 0 states, we have established the 
pair identity and the smooth variation rules applicable for 
H 2 0 lines involving high j states within individually 
defined groups such that all their spectroscopic parameters 
(i.e., the transition wavenumber, intensity, pressure broa- 
dened half-width, induced shift, and temperature expo- 
nent) must follow them [23,24]. We present our calculated 
shift values of H 2 0 lines in three groups of O'o/ <-/Tj», 
J'i/W'V}. 0/,o<— j'f.i <— and 0 V 3/-2<-J"o/'. 

j' 2 f — 2 ■*— /'i/'} in the R branch, two groups of 0/,o <-/'/' -i,i. 
j /,1 < -j"j"-i, 2 } and (| 1 2/-2^Vi,JV-2 , -/'2g’-i) in the Q. 
branch, and one group of 0 v 2 /- 2 <-j"i/' and j' 3 /- 2 <-j"o/'} 
in the P branch in Fig. 8(a)-(f), respectively. Meanwhile, 
corresponding lineshift values listed in H1TRAN 2008 are 
also presented there. According to our previous studies 
[23,24], the higher the j" is, the better the rules hold. For 
convenience, one can introduce j M , a low boundary of j", 
for a specified group as a numerical measure that the two 
rules are applicable for lines in this group whose j" values 
are above the boundary. For the six groups considered 
here, their suggested boundaries are 7, 3, 13, 5, 13, and 14, 
respectively [23,24], In Fig. 8(a)-(f), we present the calcu- 
lated shifts of lines whose/' values start from these j bd or a 
little below and up to the highest ones appearing in the list 
of the pure rotational band in H1TRAN 2008. 

As shown in Figs. 8(a) and (b) for the first two groups 
in the R branch, the calculated shifts with j" >j b d follow 
the pair identity and the smooth variation rules very well. 
On the other hand, with respect to those listed in HITRAN 
2008, values of two paired lines with j" > j bd are almost 
identical, but they could scatter widely from one pair to 
its successive pair. This implies that these values follow 
the first rule, but severely violate the second rule. We 
note that according to the reference sources in HITRAN 
2008, all of them come from the same source code 37 
denoted as recent theoretical calculations [14]. 

Then, we consider shifts of lines in the third group of 
the R branch plotted in Fig. 8(c). In general, the calculated 
shifts follow the two rules well. But, it seems they begin 
to follow the first rule with a slight delay, starting from 
j" > 14 instead of the suggested j" >13. We will explain 
why this delay could happen later. Concerning values 
provided in HITRAN 2008, some come from measure- 
ments and others from theoretical calculations. In Table 2, 
we list their reference sources provided by HITRAN 2008. 


As shown in Fig. 8(c), except for four pairs of lines with 
j"= 16, 18, 19 and 20, those values with j" >jbd in HITRAN 
2008 violate the pair identity rule. In addition, especially 
those with/' > 17 do not follow the smooth variation rule 
as well. Given the fact that usually the lineshifts are very 
small, of the order of only a few percent of the half- 
widths, to measure lineshifts are more difficult than half- 
widths. In addition, the lines with high j" are very weak, 
and to perform lineshift measurements becomes a bigger 
challenge and, therefore, measured values contain large 
uncertainties. As a result, it is understandable that mea- 
sured data do not follow the rules well. Besides, mixtures 
of different sources within the same groups also play a 
role to make the violations even worse here. It is better to 
adopt values from the same sources, as much as possible, 
for the spectroscopic parameters of lines within the same 
groups in developing databases. 

Next, for lines in the first group of the Q branch plotted 
in Fig. 8(d), the calculated shifts follow the two rules 
starting from j" >5 as expected. In this group, all the 
values listed in HITRAN 2008 come from the same source 
[14], i.e., the reference code 37. As shown in the figure, 
these values follow the pair identity rule, but somehow 
violate the smooth variation rule. More specifically, it 
seems that as/' <10, the values in HITRAN 2008 follow 
the same trend as ours, but the values with/' > 10 behave 
strangely. 

Finally, for lines in the second group of the Q branch 
and lines in the group of the P branch, results are 
presented in Fig. 8(e) and (f). As shown in these two 
figures, the calculated shifts follow the two rules. How- 
ever, they begin to follow the pair identity rule with a 
delay, i.e., roughly speaking starting from/' >16 instead 
of j" >13 and j" >14, respectively. Regarding the second 
rule, there is no such delay appearing here. With respect 
to the values listed in HITRAN 2008, in general they follow 
the first rule, but fail to follow the second one. We note 
that similar to Fig. 8(a), (b), and (d), all these values come 
from the reference source 37. 

In summary, the presented calculated shifts follow the 
two rules well. However, passing this test does not mean 
our calculated values are correct. In fact, unless one has 
made mistakes in deriving energy levels and wave func- 
tions for the H 2 0 molecule or made inconsistent errors 
somewhere else in calculations, calculated results from 
any self-consistent theories should automatically follow 
these rules. In contrast, those shift values listed in HITRAN 
2008 sometimes fail to follow the two rules, especially the 
second one. It is worth mentioning that except for several 
values in the third group of the R branch, all others come 
from the same theoretical calculations (i.e., the source 
code 37). As explained above, due to containing large 
uncertainties, measured data usually do not follow the 
rules well. However, for theoretically calculated shifts of 
lines involving high j states, their failure to follow the 
rules is fatal because usual approximations introduced in 
calculations such as simplified potential and trajectory 
models, lower order cut-offs, and so on would not cause 
the failure. Therefore, one should consider the failure as a 
clear and strong warming sign of existing inconsistent 
mistakes in formulas used or bugs in computer codes. 
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Fig. 8. Induced shifts for three groups of paired lines {/'o/ J' i/ ^j”oyl O/.o^/V.i and j' fA and {j' 3 j'- 2 <-j"oj" and j , 2/ . 2 in the R 

branch, two groups of 0 7 j'.o^/>-i,i and J/.i and {/ V 2/-2 j"-i and / 3j '-2 in the Q branch, and one group of {j v 2 j'- 2 < -/'ij" and 

j' 3 j ’ — 2 } in the P branch. They are presented in Fig. 8(a)-(f), respectively. In these plots, the calculated results are given by + and □ and they are 

connected by two solid color lines. Meanwhile, the values listed in HITRAN 2008 are plotted by x and A. 


Before ending this section, we would like to say a few 
words about the applicability of the rules. In general, 
based on the black box theory introduced in our previous 
study [23], one should expect that all the spectroscopic 
parameters of H 2 0 lines involving high j states follow 


these two rules. However, how high is high enough to 
guarantee the applicability of the two rules could differ 
slightly between the two rules themselves and among the 
different spectroscopic parameters of interest. As a con- 
venient numerical measure, the so called boundary j bd is 
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Table 2 

Reference sources of the third group in the R branch. 


Lines 

Sources* 

Lines 

Sources* 

Lines 

Sources* 

NJ 

0 

1 
O 

37 

1 63,14 *“ 1 5(3,15 

29 

203,18 *“ 1 9o,19 

37 

122 , 10 *“ 1 1 1.11 

29 

162.14*“ 15i,i5 

29 

20 2 ,i8*-19i,i9 

37 

133,11 *“ 12o,12 

29 

1 73.15 *-16o.i6 

6 

213,19 *-20o,2o 

37 

132,n<-12i,i2 

6 

172.15 *“16i,16 

6 

2l2,19*“20i,2O 

37 

14 3 .,2-13o.i3 

6 

1 83,16 *-17o,i7 

6 



14 2 .,2^13 1p13 

29 

1 82,16 *-17i,i7 

6 



153,13 <-14o,14 

29 

1 93.17 *-18o,i8 

37 



152.13*“ 14i, 14 

29 

192,17^-18, .,8 

37 




* 6: R. A. Toth, “Linelist of water vapor parameters from 500 to 8000 cm measured values, see http://mark 4 sun.jpl.nasa.gov/data/spec/H 20 . 29: R. 
A. Toth, “Air- and N 2 -Broadening parameters of water vapor: 604 to 2271 cm ! ,"J. Mol. Spectrosc. 201, 218-243 (2000). 37: Ref. [14]. 


estimated from the identity between two wave functions 
associated with paired H 2 0 states in the same category. As 
a result, there is room to adjust slightly. In addition, with 
the black box theory, one focuses attentions only on the 
relations between the inputs and outputs and does not 
address what really happens inside the box. Therefore, one 
has to consider conclusions draw from the black box theory 
to be taken in a general sense. In other words, slightly 
different behaviors to follow the general rules among 
different spectroscopic parameters are well expected. It 
turns out that spectroscopic parameters with higher 
accuracies follow the rules better than those with poorer 
accuracies. In fact, the transition positions follow these 
rules the best because this parameter is completely deter- 
mined by the identity and similarity of the energy levels of 
the H 2 0 states without involving their wave functions. 
Meanwhile, in comparison with the calculated half-widths, 
the calculated shifts begin to follow the pair identity rule 
slightly later. A theoretical explanation about this is that 
there are larger cancelations among contributions from 
different correlations to the ImS 2 (r c ) terms than those to 
ReS 2 (r c ). As a result, relative differences between ImS 2 (r c ) 
associated with two paired lines are enhanced. The above is 
a brief theoretical explanation why there are delays of the 
pair identity rule happening in Fig. 8(c), (e), and (f). 

5. Discussions and conclusions 

As discussed in our previous paper [23], all the spectro- 
scopic parameters of H 2 0 lines involving high j states, no 
matter if they are measured data or theoretically calcu- 
lated values, must follow the pair identity and smooth 
variation rules. In addition, calculated values from any 
self-consistent theories should automatically follow the 
rules. Thus, these rules can be used to check whether the 
theory used in calculations contains fatal mistakes or not. 
Our calculated half-widths follow the rules very well as 
demonstrated by our previous works [23,24] and mean- 
while, as shown in Fig. 8, our calculated shifts follow the 
rules also. In contrast, some of calculated shift values in 
HITRAN 2008 fail to follow the rules. A question arises that 
how this could happen and where the mistakes could be? 

We would like to note that this question differs from 
answering why usually the accuracy qualities of calculated 
shifts are poorer than those of calculated half-widths. The 
former is searching for mistakes made by people in their 
calculations and the latter is explaining different natures 


between the half-width and the lineshift. However, the 
ways to solve the problems are the same: one has to 
investigate the entire process within the specified formal- 
ism to figure out the answers. Because the formalism used 
in our calculations differ from the one used in calculating 
theoretical results in HITRAN 2008 [14] and because we 
have no access to their codes, we do not know exactly the 
answers at this point. But, based on our experiences, we do 
know where to search for them and we would like to 
provide some suggestions here. The first place required to 
check is the energy levels and wave functions of H 2 0 
states, especially those associated with high j values. If 
calculated energy levels are not accurate enough, one 
should replace them by more accurate ones available in 
literature [25]. With respect to the wave functions, one 
should pay attention to their symmetries associated with 
the four sub-blocks (i.e., E + , E _ , 0 + , and O - ) [23,24] 
because violations of these symmetries are not allowed. In 
addition, we would like to note that by selecting the 
representation alone, the orientation of H 2 0 in the mole- 
cular-fixed frame has not been uniquely determined. For 
example, in the I R representation where the principle 
inertial axis b of H 2 0 is coincided with the x axis of the 
molecular fixed frame, the dipole moment could still be 
either positive or negative because the O atom can be 
located at either the negative or the positive sides of the x 
axis. Therefore, one has to make sure the choices of the O 
atom location are the same both in developing the wave 
functions and in describing the interaction potentials of the 
H 2 0-N 2 pair. Given the fact that subroutines used in 
calculating the wave functions may come from other 
people, it becomes necessary to check these to avoid any 
inconsistency. 

Meanwhile, one should check the resonance functions 
carefully. As long as potential models contain the short- 
range atom-atom component, one has to evaluate two 
sets of resonance functions consisting of huge members. 
The one set used in calculating R eS 2 (r c ) consists of 
resonance functions whose magnitude decrease very 
quickly as their arguments increase. In contrast, the other 
set used in calculating ImS 2 (r c ) consists of resonance 
functions whose magnitudes decrease very slowly as their 
arguments go to infinity. This implies that one has to 
evaluate them over wider ranges. As a result, to evaluate 
the second set of the resonance functions is a big 
challenge. We suspect that not enough attention has been 
paid to the different behaviors between these two sets of 
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resonance functions and that the second set has not been 
well depicted over their whole ranges which could intro- 
duce significant errors for calculated shifts. 

Now, we summarize the salient features of the present 
study. By adopting the concept of the causal function from 
signal processing and relying on the sampling theory to 
replace the continuous Fourier transform by the discrete 
Fourier transform, we have discovered powerful tools 
which work excellently in evaluating the Hilbert trans- 
forms without performing the Cauchy principal integra- 
tions. With this new method, we are able to effectively 
and accurately calculate converged values of the N 2 
induced shifts of H 2 0 lines. Thus, the challenge to calcu- 
late converged lineshifts with our formalism developed 
using the coordinate representation has finally been over- 
come. As a result, one is able to calculate both pressure 
broadened half-widths and pressure induced shifts to the 
accuracy of the approximations in the interaction-poten- 
tial and trajectory models without containing conver- 
gence errors within the current frameworks of the 
original and the modified RB formalism. 

Furthermore, after we were able to derive the Hilbert 
transforms, we have achieved fruitful results by finding 
answers to several longstanding questions. By comparing 
our results with those listed in HITRAN 2008, most of 
which are theoretically calculated values using the same 
potential model, one can judge how large their differences 
are. Then, by screening both calculated results with the 
two rules, one can clearly identify which ones contain fatal 
errors. Secondly, we are able to demonstrate significant 
effects on calculated shifts from the modification of the 
original RB formalism. Finally, we are able to verify the 
validity of the assumption with which one calculates the 
half-widths without including contributions from lmS 2 (r c ). 

The next step necessary to improve theoretically calcu- 
lated results is using a properly selected set of the most 
accurately measured experimental values of the half-width 
to refine the potential parameters and then calculate 
calculated half-width values for comparison. Based on 
the new optimized potential model, one can then calculate 
theoretical values of half-widths and shifts for spectral 
lines of interest. Meanwhile, one has to keep in mind that 
some fundamental problems are required to address in the 
current RB formalism, such as how to consider the cou- 
plings between the translational and internal motions and 
so on. Therefore, there are new challenges ahead. 
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